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Recent analysis of the isospin diffusion data from heavy- ion collisions based on an isospin- 
and momentum-dependent transport model with in-medium nucleon-nucleon cross sec- 
tions has led to the extraction of a value of L = 88 ± 25 MeV for the slope of the nuclear 
symmetry energy at saturation density. This imposes stringent constraints on both the 
parameters in the Skyrme effective interactions and the neutron skin thickness of heavy 
nuclei. Among the 21 sets of Skyrme interactions commonly used in nuclear structure 
studies, the 4 sets SIV, SV, G CT , and R CT are found to give L values that are consistent 
with the extracted one. Further study on the correlations between the thickness of the 
neutron skin in finite nuclei and the nuclear matter symmetry energy in the Skyrme 
Hartree-Fock approach leads to predicted thickness of the neutron skin of 0.22 ± 0.04 fm 
for 208 Pb, 0.29 ± 0.04 fm for 132 Sn, and 0.22 ± 0.04 fm for 124 Sn. 



1. Introduction 

The study of the equation of state (EOS) of isospin asymmetric nuclear matter, 
especially the nuclear symmetry energy, is currently an active field of research 
in nuclear physics i- 2 - 3 ' 4 ' 5 . 6 ' 7 ' 8 . Although the nuclear symmetry energy at nor- 
mal nuclear matter density is known to be around 30 MeV from the empirical 
liquid-drop mass formula 9,1 °, its values at other densities, especially at supra- 
normal densities, are poorly known ' . Advances in radioactive nuclear beam 
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facilities provide, however, the possibility to pin down the density dependence 
of the nuclear symmetry energy in heavy ion collisions induced by these nuclei 

1,2,7,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35 J nc J eec J significant 

progress has recently been made in extracting the information on the density depen- 
dence of nuclear symmetry energy from the isospin diffusion data in heavy-ion colli- 
sions at NSCL/MSU 36 < 37 < 38 . Using an isospin- and momentum-dependent IBUU04 
transport model with in-medium nucleon-nucleon (NN) cross sections, the isospin 
diffusion data were found to be consistent with a relatively soft nuclear symmetry 
energy at subnormal density 38 . 

Information on the density dependence of the nuclear symmetry energy can also 
be obtained from the thickness of the neutron skin in heavy nuclei as the latter is 
strongly correlated with the slope L of the nuclear matter symmetry energy at 
saturation density 39, 40,41,42, 43, 4_ B ecause f the large uncertainties in measured 
neutron skin thickness of heavy nuclei, this has not been possible. Instead, studies 
have been carried out to use the extracted nuclear symmetry energy from the isospin 
diffusion data to constrain the neutron skin thickness of heavy nuclei 44 ' 38 . Using 
the Hartree-Fock approximation with parameters fitted to the phenomenological 
EOS that was used in the IBUU04 transport model to describe the isospin diffusion 
data from NSCL/MSU, it was found that a neutron skin thickness of less than 0.15 
fm 44 ' 38 for 208 Pb was incompatible with the isospin diffusion data. 

In the present talk, we report our recent work on constraining the Skyrmc 
effective interactions and the neutron skin thickness of heavy nuclei using the isospin 
diffusion data from heavy ion collisions 45 . Using the value of L obtained from 
the extracted density dependence of the nuclear symmetry energy from the isospin 
diffusion data, we have been able to limit the allowed parameter sets for the Skyrme 
interaction. Also, studying the correlation between the density dependence of the 
nuclear symmetry energy and the thickness of the neutron skin in a number of 
nuclei within the framework of the Skyrmc Hartree-Fock approach further allows 
us to obtain stringent constraints on the neutron skin thickness of the nuclei 208 Pb, 
132 Sn, and 124 Sn. 



2. Nuclear symmetry energy and the Skyrme interaction 

The nuclear symmetry energy E sym (p) at nuclear density p can be expanded around 
the nuclear matter saturation density p n as 



where L and K synl are the slope and curvature of the nuclear symmetry energy at 
Pa, i-e., 

T _„ dE sym (p) _ 2 d 2 E sym {p) 

~df) l' 9=p0 ' A sym — y/?o \P=Po- \ l ) 
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The L and K syra characterize the density dependence of the nuclear symmetry en- 
ergy around normal nuclear matter density, and thus provide important information 
on the properties of nuclear symmetry energy at both high and low densities. 

In the standard Skyrme Hartree-Fock approach 46 > 4 7,48,49,50,39^ the interaction is 
taken to be zero-range and density- and momentum-dependent with its parameters 
fitted to the binding energies and charge radii of a large number of nuclei in the 
periodic table. For infinite nuclear matter, the nuclear symmetry energy from the 
Skyrme interaction can be expressed as 49 > 50 

Esym(p) = 3 ^ [-) P V3 - ^t (2x + l)p - —t 3 (2x 3 + 

+ Ya (¥) ' hMlXl + {4 + 5X2) <2]p5/3 ' (3) 

where the a, to — 13, and xq — x% are Skyrme interaction parameters. 




Fig. 1. (Color online) Density dependence of the nuclear symmetry energy -E sym (p) for 21 
sets of Skyrme interaction parameters (solid lines) as well as the MDI interaction with 
x = — 1 (open squares) and (solid squares). 



Fig. 1 displays the density dependence of E sym (p) for 21 sets of Skyrme interac- 
tion parameters, i.e., SKM, SKM*, RATP, SI, SII, SIII, SIV, SV, SVI, E, E CT , G CT , 
R CT , Z, Zo-, Z*, T, T3, SkX, SkXce, and SkXm. Values of the parameters in these 
Skyrme interactions can be found in Refs. 46 > 47 > 48 . For comparison, we also show 
in Fig. 1 results from the phenomenological parametrization of the momentum- 
dependent nuclear mean-field potential based on the Gogny effective interaction 51 , 
i.e., the MDI interactions with x = — 1 (open squares) and (solid squares), where 
different x values correspond to different density dependence of the nuclear symme- 
try energy but keep other properties of the nuclear EOS the same . By comparing 
the isospin diffusion data from NSCL/MSU with results from the IBUU04 trans- 
port model using in-medium NN cross sections, these interactions have been shown 
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to give, respectively, the upper and lower bounds for the stiffness of the nuclear 
symmetry energy 38 . It is seen from Fig. 1 that the density dependence of the sym- 
metry energy varies drastically among different interactions. Although the values 
of -Esym(po) are all in the small range of 26-35 MeV, the values of L and K aym are 
in a much larger range of —50-100 MeV and —700-50 MeV, respectively. 

3. Constraining symmetry energy from isospin diffusion data in 
heavy ion collisions 

Nuclear symmetry energy is known to affect the diffusion of isospin in intermediate- 
energy heavy ion collisions. Experimentally, the degree of isospin diffusion between 
the projectile nucleus A and the target nucleus B can be studied via the quantity 

n 52,36 

2X A+B _ X A+A _ X B+B 

* ~ X A + A - x B+B ' 
where X is any isospin-sensitive observable. By construction, the value of Ri is 
1 (—1) for symmetric A + A (B + B) reaction. If isospin equilibrium is reached 
during the collision as a result of isospin diffusion, the value of Ri is about zero. In 
the NSCL/MSU experiments with A = 124 Sn and B = 112 Sn at a beam energy of 
50 MeV/nucleon and an impact parameter about 6 fm, the isospin asymmetry of 
the projectile-like residue was used as the isospin tracer X 36 . Using an isospin- and 
momentum-dependent IBUU04 transport model with free-space experimental NN 
cross sections, the dependence of Ri on the nuclear symmetry energy was studied 
from the average isospin asymmetry of the projectile-like residue that was calculated 
from nucleons with local densities higher than po/20 and velocities larger than 1/2 
the beam velocity in the center-of-mass frame 37 . Comparing theoretical results 
with experimental data has allowed us to extract a nuclear symmetry energy of 
E sym (p) ~ 31.6(p/po) 105 - Including also medium-dependent NN cross sections, 
which are important for isospin-dependent observables 38 > 53 J the isospin diffusion 
data leads to an even softer nuclear symmetry energy of E sym (p) w 31.6(p/ po) 7 
with 7 w 0.7 38 . 

In Fig. 2, we show the degree of the isospin diffusion 1 — Ri as a function of L 
obtained from the IBUU04 transport model with in-medium NN cross sections and 
the mean-field potential based on the MDI interactions. The shaded band in Fig. 2 
indicates the data from NSCL/MSU 36 . It is seen that the strength of isospin diffu- 
sion 1 — Ri decreases monotonically with decreasing value of x or increasing value 
of L. This is expected as the parameter L reflects the difference in the pressures on 
neutrons and protons. From comparison of the theoretical results with the data, we 
can clearly exclude the MDI interaction with x = 1 and x — — 2 as they give either 
too large or too small a value for 1 — Ri compared to that of data. The range of x 
or L values that give values of 1 — Ri falling within the band of experimental values 
could in principle be determined in our model by detailed calculations. Instead, 
we determine this schematically by using the results from the four x values. For 
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the centroid value of L, it is obtained from the interception of the line connecting 
the theoretical results at x = — 1 and with the central value of 1 — Ri data in 
Fig. 2, i.e., L = 88 MeV. The upper limit of L = 113 MeV is estimated from the 
interception of the line connecting the upper error bars of the theoretical results at 
x = — 1 and —2 with the lower limit of the data band of 1 — R4. Similarly, the lower 
limit of L = 65 MeV is estimated from the interception of the line connecting the 
lower error bars of the theoretical results at x — and —1 with the upper limit of 
the data band of 1 — B4. This leads to an extracted value of L = 88 ± 25 MeV as 
shown by the solid square with error bar in Fig. 2. 

1.0 

0.8 

Of 0.6 
0.4 



0.2 

30 60 90 120 150 180 
L (MeV) 

Fig. 2. (Color online) Degree of the isospin diffusion 1 — Ri as a function of L using the MDI 
interaction with x = —2, —1, 0, and 1. The shaded band indicates data from NSCL/MSU 
36 . The solid square with error bar represents L = 88 ± 25 MeV. 

The extracted value of L — 88±25 MeV gives a rather stringent constraint on the 
density dependence of the nuclear symmetry energy and thus puts strong limitations 
on the nuclear effective interactions as well. For the Skyrme effective interactions 
shown in Fig. 1, for instance, all of those lie beyond x — and x = — 1 in the 
sub-saturation region are not consistent with the extracted value of L. Actually, we 
note that only 4 sets of Skyrme interactions, i.e., SIV, SV, G CT , and R CT , in the 21 
sets of Skyrme interactions considered here have nuclear symmetry energies that 
are consistent with the extracted L value. 

4. Correlations between neutron skin thickness of finite nuclei and 
the nuclear symmetry energy at saturation density 

Also affected by the density dependence of nuclear symmetry energy is the neutron 
skin thickness S of a nucleus, which is defined as the difference between the root- 
mean-square radii yj (r n ) of neutrons and \J (r p ) of protons, i.e., 

s = VWi-J(jf)- (5) 



x= 1 

x=0 


E/A=50 MeV and b=B fm 
-□- IBUU04 with MDI + a ma 

L=88t25 MeV 

v 
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In particular, S is sensitive to the slope parameter L of the nuclear symmetry energy 
at normal nuclear matter density 39 > 40 > 41 > 42 > 43 > 4 . Using above 21 sets of Skyrme 
interaction parameters, we have evaluated the neutron skin thickness of several 
nuclei. In Figs. 3(a), (b) and (c), we show, respectively, the correlations between 
the neutron skin thickness of 208 Pb with L, K synl , and E sym (p ). It is seen from 
Fig. 3(a) that there exists an approximate linear correlation between S and L. The 
correlations of S with K sym and E sym (p ) are less strong and even exhibit some 
irregular behavior. The solid line in Fig. 3(a) is a linear fit to the correlation between 
S and L and is given by the following expression: 

S( 208 Pb) = (0.1066 ± 0.0019) + (0.00133 ± 3.76 x 10~ 5 ) x L, (6) 

or 

L = (-78.5 ± 3.2) + (740.4 ± 20.9) x S^Pb), (7) 

where the units of L and S are MeV and fm, respectively. Therefore, if the value 
for either S'( 208 Pb) or L is known, the value for the other can be determined. 
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Fig. 3. (Color online) Neutron skin thickness S of Pb as a function of (a) L, (b) -Ksym, 
and (c) -Esym(po) f° r 21 sets of Skyrme interaction parameters. The line in panel (a) 
represents a linear fit. 



It is of interest to see if there are also correlations between the neutron skin 
thickness of other neutron-rich nuclei and the nuclear symmetry energy. Fig. 4 
shows the same correlations as in Fig. 3 but for the neutron-rich nuclei 132 Sn, 
124 Sn, and 48 Ca. For the heavier 132 Sn and 124 Sn, we obtain a similar conclusion as 
for 208 Pb, namely, S exhibits an approximate linear correlation with L but weaker 
correlations with K sym and E sym (po). For the lighter 48 Ca, on the other hand, 
all the correlations become weaker than those of heavier nuclei. The neutron skin 
thickness of heavy nuclei is thus better correlated with the density dependence of 
the nuclear symmetry energy. As in Eq. (6) and (7), a linear fit to the correlation 
between S and L can also be obtained for 132 Sn and 124 Sn, and the corresponding 
expressions are 

S( 132 Sn) = (0.1694 ± 0.0025) + (0.0014 ± 5.12 x 10~ 5 ) x L, (8) 
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Table 1. Linear correlation coefficients C; of S with 
L, K sym and £ sym (po) for 208 Pb, 132 Sn, 124 Sn, and 
48 Ca from 21 sets of Skyrme interaction parameters. 



c, (%) 


208p b 


132 Sn 


124 Sn 


48 Ca 


S-L 


99.25 


98.76 


98.75 


93.66 


S-Ksym 


92.26 


92.06 


92.22 


86.99 


S-E S y m 


87.89 


85.74 


85.77 


81.01 



L = (-117.1 ± 5.4) + (695.1 ± 25.3) x 5( 132 Sn), (9) 

and 

S*( 124 Sn) = (0.1255 ± 0.0020) + (0.0011 ± 4.05 x 10~ 5 ) x L, (10) 

L = (-110.1 ± 5.2) + (882.6 ± 32.3) x S( 124 Sil), (11) 
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Fig. 4. (Color online) Same as Fig. 3 but for nuclei Sn (Solid squares), Sn (Open 
squares) and 48 Ca (Triangles). 



Similar linear relations between S and L are also expected for other heavy nuclei. 
This is not surprising as detailed discussions in Refs. 39 > 4 o, 41,42,43, 4 have shown that 
the thickness of the neutron skin in heavy nuclei is determined by the pressure 
difference between neutrons and protons, which is proportional to the parameter 
L. 

To give a quantitative estimate of above discussed correlations, we define the 
following linear correlation coefficient Cf. 

d = y/l - q/t, (12) 

where 

n n n 

i— 1 i— 1 i— 1 

In the above, A and B are the linear regression coefficients, {pa, yi) are the sample 
points, and n is the number of sample points. The linear correlation coefficient C; 
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measures the degree of linear correlation, and Ci — 1 corresponds to an ideal linear 
correlation. Table 1 gives the linear correlation coefficient C; for the correlation of 
S with L, K syra and E sym (p a ) for 208 Pb, 132 Sn, 124 Sn, and 48 Ca shown in Figs. 3 
and 4 for different Skyrmc interactions. It is seen that these correlations become 
weaker with decreasing nucleus mass, and a strong linear correlation only exists 
between the S and L for the heavier nuclei 208 Pb, 132 Sn, and 124 Sn. Therefore, the 
neutron skin thickness of these nuclei can be extracted once the slope parameter L 
of the nuclear symmetry energy at saturation density is known. 

5. Predictions on the neutron skin thickness of heavy nuclei 

The extracted L value from the isospin diffusion data in heavy ion collisions allows 
us to determine from Eqs. (6), (8), and (10), respectively, a neutron skin thickness 
of 0.22 ± 0.04 fm for 208 Pb, 0.29 ± 0.04 fm for 132 Sn, and 0.22 ± 0.04 fm for 124 Sn. 
Experimentally, great efforts were devoted to measure the thickness of the neutron 
skin in heavy nuclei 54 > 55 , and a recent review can be found in Ref. 56 . The data 
for the neutron skin thickness of 208 Pb indicate a large uncertainty, i.e., 0.1-0.28 
fm. Our results for the neutron skin thickness of 208 Pb are thus consistent with 
present data but give a much stronger constraint. A large uncertainty is also found 
experimentally in the neutron skin thickness of 124 Sn, i.e., its value varies from 0.1 
fm to 0.3 fm depending on the experimental method. The proposed experiment of 
parity-violating electron scattering from 208 Pb at the Jefferson Laboratory is ex- 
pected to give another independent and more accurate measurement of its neutron 
skin thickness (within 0.05 fm), thus providing improved constraints on the density 
dependence of the nuclear symmetry energy 57 ' 58 . 

Most recently, an accurately calibrated relativistic parametrization based on 
the relativistic mean-field theory has been introduced to study the neutron skin 
thickness of finite nuclei 59 . This parametrization can describe simultaneously the 
ground state properties of finite nuclei and their monopole and dipole resonances. 
Using this parametrization, the authors predicted a neutron skin thickness of 0.21 
fm in 208 Pb, 0.27 fm in 132 Sn, and 0.19 fm in 124 Sn 59 < 60 . These predictions are in 
surprisingly good agreement with our results constrained by the isospin diffusion 
data from heavy-ion collisions. 

6. Summary 

Both the structure of nuclear surface and the dynamics of isospin equilibration 
in heavy ion collisions are affected by the density dependence of nuclear symmetry 
energy. From the most recent analysis of the isospin diffusion data in heavy-ion colli- 
sions using an isospin- and momentum-dependent transport model with in-medium 
NN cross sections, a value L = 88 ± 25 MeV has been extracted for the slope of the 
nuclear symmetry energy at saturation density. This relatively constrained value 
imposes strong constraints on the parameters in the Skyrme effective interactions. 
Among the 21 sets of commonly used Skyrme parameters, only SIV, SV, G a , and 
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R a have L values that are consistent with the extracted one. We have also studied 
the correlation between the neutron skin thickness of finite nuclei and the nuclear 
symmetry energy within the framework of the Skyrme Hartree-Fock model. As in 
previous studies, we have found a strong linear correlation between the neutron skin 
thickness of heavy nuclei and the slope L of the nuclear matter symmetry energy at 
saturation density. This correlation provides stringent constraints on both the den- 
sity dependence of the nuclear symmetry energy and the thickness of the neutron 
skin in heavy nuclei. The extracted L value from the isospin diffusion data then 
leads to predicted neutron skin thickness of 0.22 ± 0.04 fm for 208 Pb, 0.29 ± 0.04 
fm for 132 Sn, and 0.22 ± 0.04 fm for 124 Sn. Our work has thus demonstrated that 
information obtained from studying isospin dynamics in heavy ion collisions can be 
very useful for understanding the structure of nuclei. 
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